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A Method for transforming a digital signal from the time 

domain to the frequency domain 

Background of the Invention 

5 

Domain transformation, specifically Discrete cosine transforms (DCT), are widely used 
in modern signal processing industry. Recently, a 'relative' of the DCT, called integer 
DCT, has attracted a lot of research interests because of its important role in lossless 
coding applications. The term 'lossless' means that the decoder can generate an exact 
10 copy of the source signal from the encoded bit-stream. 

The DCT is a real-valued block transform. Even if the input block consists of all integers, 
the output block of the DCT are made of real numbers. For convenience, the input block 
is referred as input vector and the output block as output vector. If a vector has all integer 
15 components, it is called an integer vector. In contrast to the DCT, the integer DCT 
generates an integer output vector from an integer input vector. For the same integer 
input vector, the integer output vector of integer DCT closely approximates the real 
output vector of DCT. Thus the integer DCT keeps all the good properties of the DCT in 
spectrum analysis. 

20 

An important property of the integer DCT is reversibility. Reversibility means that there 
exists an integer inverse DCT (IDCT) so that if the integer DCT generates output vector y 
from input vector jc, the integer IDCT can recover vector jc from vector^. Sometimes the 
integer DCT is also referred as the forward transform, and the integer IDCT as the 

25 backward or inverse transform. 

i 

A transform called integer modified discrete cosine transform (IntMDCT) is recently 
proposed and used in the ISO/IEC MPEG-4 audio compression. It is shown in [2] that the 
core of IntMDCT is an integer type-IV DCT (DCT-IV). This motivates us to search for 
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an efficient fast algorithm, specifically for the integer DCT-IV algorithm. There are four 
types of DCT, and the fourth type is of particular interest. Fast algorithms refer to those 
efficient methods to compute a transform. 

To obtain a fast algorithm for integer DCT-IV, the direct approach is to modify existing 
DCT-IV fast algorithms. As an example, Wang gave a DCT-IV fast algorithm in [4]. This 
algorithm can be converted to an integer DCT-IV fast algorithm by factorizing each 
butterfly operation (Givens rotation) into three 'lifting steps*. A detailed explanation of 
this factorization method can be found in [2]. For this type of integer DCT-IV algorithms, 
the total rounding number is at level N\og 2 N , where N is the block size. Here, rounding 
means the operation to round a real number to its nearest integer. Due to roundings, there 
is a difference between the outputs of integer DCT-IV and DCT-IV under same inputs. 
This difference is referred as approximation error. Generally, this approximation error 
increases with the total rounding number of the integer algorithm. A good integer DCT- 
IV algorithm should use as less roundings as possible. 

r 

In this invention, a method for transforming a digital signal from the time domain to the 
frequency domain, in particular reversible integer transform like DCT IV is introduced. 
The total rounding number of the method according to the invention can be significantly 
reduced, for example in case of DCT IV to be as low as 1.5N. As a result, the 
approximation error of the method according to the invention is far less than that of the 
directly converted integer transforms. The computational complexity of die proposed 
method is also very low. 

Summary of the Invention ^ 

The invention relates to a method for transforming a digital signal from the time domain 
to the frequency domain. The method can be used for any types of digital signal, such as 
audio, image or video signals. The digital signal, which corresponds to a physical 
measured signal, may be generated by scanning at least a characteristic feature of a 
corresponding-analog signal (for example, the luminance and chrominance values of a 
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video signal, the amplitude of an analog sound signal, or the analog sensing signal from a 
sensor). The digital signal comprises a plurality of data symbols. The data symbols of 
the digital signal are grouped into blocks, with each block having the same predefined 
number of data symbols based on the sampling rate of the corresponding analog signal. 

5 

The digital signal is transformed to the frequency domain by a transforming element 
based on a transformation function which is implemented within a transformation matrix. 
In this method according to the invention, the transforming element transforms the digital 
signal by processing two blocks of data symbols as input signals simultaneously to 
10 generate two corresponding blocks of transformed output signal. It should be noted that 
the block size of the input signals and the output signals are the same. 

An example of a transformation matrix is: 



15 



c w 



wherein C N is a matrix representing the transformation function. shall be referred 
as the transformation function matrix henceforth. 



20 The transforming element includes a plurality of lifting stages. A first lifting stage 
receives the two sub-blocks as first and second input signals, respectively. The first input 
signal is processed in a transformation path in the first lifting stage which includes a 
domain transformer, a rounding unit and a summation unit. The domain transformer can 
be any type of non-integer (real value) transformation function. The transformed first 

25 input signal is summed with the second input signal and result in an output signal. The 
output signal and the first input signal of the first lifting stage are used as the first and 
second input signals for a second (or subsequent) lifting stage, resulting in an output 
signal by the second lifting stage. Similarly, the output signal and the first input signal of 
the second lifting stage are used as the first and second input signals for another 

30 subsequent lifting stage. Since the output signal and the first input signal from each stage 
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are received by the subsequent stage, these two signals shall be considered to be signals 
output from each of the lifting stages, although the second input signal of the subsequent 
stage is identical to the first input signal of die previous stage (i.e. not transformed). 

5 The transforming element can be illustrated based on the model of a lifting ladder. The 
lifting ladder model has two side pieces, each for receiving one of the two blocks of data 
symbols. Two or more cascading lifting stages are provided between the two side pieces. 
Each lifting stage receives a signal at one end (input end), and outputs a signal at the 
other end (output end) via a summation unit A domain transformer is arranged at the 
10 input end of the lifting stage, and a rounding unit is arranged at the output end, between 
the domain transformer and the summation unit. The lifting stages are arranged between 
the side pieces in an alternating manner, such that the output (or input) ends of adjacent 
lifting stages are connected to the different side pieces. 

« 

15 It should be noted that although the transforming element is described in the form of the 
lifting ladder model, it is only to illustrate the transformation paths of the transforming 
element. However, the invention shall not be limited to said ladder model. 

Figure 1 shows a flow chart of an embodiment of the method according to the invention 
20 using three lifting stages. In Fig.l, xj and X2 are first and second blocks of the digital 
signal, respectively, z is an intermediate signal, and yj and yi are output signals 
corresponding to the respective first and second block of the digital signal. 

The method according to the invention can be used for transforming an input digital 
25 signal which represents integer values to an output signal which also represents integer 
values. In other words, a "same word-length" transformation is performed such that 
word-lengths of the input signal and the output signal are the same (for example, a 8-bit 
input data is transformed to a 8-bit output data). Since both the input and output signals 
have the same word-length, the transformation method according to the invention is 
3D xeversible. The output signal may be transformed back to the original input signal by 
performing the transformation method according to the invention. Such a reversibility 
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property of the transformation according to the method of the invention can be used in 
lossless coding in which the output signal should be identical to the original input signal. 

Due to this reversibility property, the invention shall also include the transforming of a 
5 digital signal from the frequency domain to the time domain using the transformation 
element as described earlier. 

Such a same word-length (or integer) transformation of signals according to the invention 
can be used in many applications and systems such as MPEG audio, image and video 
10 compression, JPEG2000 or spectral analysers (for analyzing Infrared, Ultra-violet or 
Nuclear Magnetic Radiation signals). It can also be easily implemented in hardware 
systems such as in a fixed-point Digital Signal Processor (DSP), without having to 
consider factors such as overflow in the case of a real-value signal transformation. 

15 The transformation function used for the transformation of the digital signal according to 
the invention can be any kind of transformation function which can be represented as an 
Involutory Matrix. An Involutory Matrix is a matrix when multiplied with itself forms an 
Identity Matrix. 

20 Let Tbe a matrix. T is an Involutory Matrix if it satisfies the following equation: 
T'T = I„, 



25 



where/ 



N 



1 



1 



is an Identity Matrix. 



i 

Examples of transformation functions which can be represented as an Involutory Matrix 

» 

includes, but not limited to, Discrete Cosine Transform (DCT) I, DCT IV, Discrete Sine 
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Transform (DST) I, DST IV, Discrete Fourier Transform (DFT) I, DFT IV, Discrete 
Wavelet Transform (DWT) I and Discrete Wavelet Transform IV. Therefore, these 
transformation functions can be used as the transformation function in the method 
according to the invention. 

5 

The domain transformers in the transforming element are implemented using the same 
transformation function that are used within the transformation matrix. However, the fast 
algorithms for implementing the transformation function for each of the domain 
transformers may be different. For example, when DCT IV is used as the transformation 

10 function for transforming the input digital signal, all the domain transformers should also 
be implemented using a DCT IV transformation function. However, various algorithms 
can be used to implement the DCT IV transformation function. Such algorithms include 
DCT II or FFT algorithms. Therefore, a different algorithm may be used in each of the 
domain transformers for implementing the DCT IV transformation functions. 

15 Alternatively, the same algorithm may be used in all the domain transformers. 

As can be seen from the method according to the invention, all the data symbols in each 
block of the digital signal are provided to the transforming element as a data vector. In 
each lifting stage, the data vector is transformed in the domain transformer, and the 

20 transformed data vector is rounded subsequently to an integer vector (i.e. after the 
transformation in each lifting stage). In other words, rounding is performed in the method 
of the invention once on the data vector as a whole. This is in contrast to any method 
according to the state of the art, wherein the rounding process is performed within the 
transformation process for the individual element or data symbol in each block of the 

25 digital signal. Thus, the number of rounding operations in the method according to the 
invention is greatly reduced. Due to the reduced number of rounding operations, the 
method according to the invention does not require large computation time and computer 
resource. 

* 

30 The transformation matrix, which comprises the transformation function as sub-matrices, 
can be further decomposed into a plurality of lifting matrices. Each of the lifting matrices 
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corresponds to each lifting stage of the transformation element. A mathematical example 
illustrating this implementation will be shown later. . 

Each lifting matrix comprises four sub-matrices, with two permutation matrices as two of 
the sub-matrices in one diagonal, and a zero and a transformation function matrix as the 
other two of the sub-matrices in the other diagonal. A permutation matrix is a matrix 
which changes the position of the elements in another matrix. An example of a lifting 
matrix is: 



10 



15 



P2 
C 



IV 



0 

PI 



wherein C% is the transformation function matrix, and P\ N and P2 N are permutation 
matrices. The permutation matrices may also be identical. The permutation matrices 



may also be arranged in the lifting matrix as 



r 0 
P2 



PI 



N 



AT 



Also, Identity Matrix is commonly used as a form of permutation matrix. An Identity 
Matrix is a matrix shown in the following: 



Is = 



1 



1 



The number of lifting matrices, and hence the number of lifting stages in the 
20 transformation element, is three, if the DCT-IV is the transformation function. 

It should be noted that each of the elements in the lifting matrices is a matrix itself. 
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The invention not only relates to a method for transforming a digital signal from the time 
domain to the frequency domain (and vice versa), but also includes a computer program, 



5 Example of the method according to the invention based on an Integer-type 
DCT IV as the transformation function. 

An example of the method according to the invention is illustrated using the DCT-IV 
transformation function. The DCT-IV of a JV-point real input sequence x(n) is defined 

10 as follows [3]: 





(1) 



Let Cjj be the transformation matrix of DCT-IV, that is 



15 




(2) 



The following relation holds for the inverse DCT-IV matrix [3]: 



20 




(3) 



Let X = [x(n)] 



n-OJ.-s/V-l 



and y = [y(m)] 



fli-»o.u—jv-i 



be two N x 1 column vectors, Equation 



(l)-can be expressed as 




(4) 
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Now, assume that there are two real, integer Nxl column vectors x v x 2 . The column 
vectors x,,x 2 correspond to the first and second blocks of the digital signal The DCT-IV 
transforms of x v x 2 are y x ,y 2 , respectively. 



5 Vl = C»x x 

y 2 =» c*n x 2 



(5) 
(6) 



Combining (5) and (6): 



10 



Y\ 

lV2 



JL^2j 



(7) 



15 



Or, alternatively, 



Y\ 

Y21 



,IV 
N 



1*1 



Let T 2N be the counter diagonal matrix of C„ in (8), that is 



(8) 



IN 



C IV1 



ilV 



(9) 



20 Matrix T 2N can be factorized as follows: 



IN 









r IV 

N 







(10) 



where 1 N is the N x N identity matrix. 



25 
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I 



Equation (10) can be easily verified using the DCT-IV property in (3). Using (10), 
Equation (8) can be expressed as 



r / 



N 



-C" I 



N *N 



J . 



N 



-/„ c 



Vat 1 n\\j*\ 



(ii) 



From (1 1), the following integer DCT-IV algorithms that simultaneously compute two 
integer DCT-IVs is derived. 

The time to frequency domain integer transform is determined by the following: 



10 



Step 1) Compute 



(12a) 



15 Step 2) Compute 



y x -[c%z\-x. 



(12b) 



20 



Step 3) Compute 



(12c) 



25 



where |*J denotes rounding operation. Figure 2 illustrates the method of transforming a 

digital signal from the time domain to the frequency domain using DCT-IV as the 
transformation function, xj and xi are two blocks of the input digital signal, z is an 
intermediate signal, eaidyj and j>2 are corresponding blocks of the output signal. 



"The frequency to time domain integer transform is determined by the following: 
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Step 1) Compute 



(13a) 



5 Step 2) Compute 



(13b) 



Step 3) Compute 



10 



+ z 



(13c) 



Figure 3 illustrates the method of transforming a digital signal from the frequency 
domain to the time domain using DCT-IV as the transformation function. 

15 

Conclusion of the above example 

Equation (12) and (13) show that to compute two N x N integer DCT-IVs, three N x N 
DCT-IVs, three Nxl roundings, and three JVxl additions are needed. Therefore, for 
20 one N x N integer DCT-IV, the average is: 

RC(N)-l.5N (14) 
AC(N) = l.5AC(C'J) + l.5N 1 (15) 



where RC( ) is the total rounding number, and AC( ) is the total number of arithmetic 
operations. Compared to the directly converted integer DCT-IV algorithms, the proposed 
integer DCT-IV algorithm reduces RC from level N\o% 2 N to N. This is the result of 

doing rounding after the DCT-IV transform instead of in the DCT-IV transform. 



25 
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As indicated by (15), the arithmetic complexity of the proposed integer DCT-IV 
algorithm is about 50 percent more man that of a DCT-IV algorithm. However, if RC is 
also considered, the combined complexity (AC+RC) of the proposed algorithm does not 
5 much exceed that of the directly converted integer algorithms. Exact analysis of the 
algorithm complexity depends on the DCT-IV algorithm used. 

As shown in Figure 1 and 2, the proposed IntDCT-IV algorithm is simple and modular in 
structure. It can use any existing DCT-IV algorithms in its DCT-IV computation block. 
10 The proposed algorithm is suitable for applications that require IntMDCT, e.g. in the 
MPEG-4 audio extension 3 reference model 0. 

» 

The proposed algorithm is suitable for both mono and stereo applications. In mono 
applications, two consecutive blocks of samples are grouped and processed together. This 
15 introduces signal delay of one block length when compared to single-block processing. 
However, in stereo applications, this extra block delay can be prevented, if simultaneous 
sample blocks from the left and the right channel are grouped and processed together. 

In this invention, a method for realizing reversible transformation function, for example 
20 for integer type-IV DCT transformation function, is proposed. This method requires 
significantly reduced number of roundings for every block of N input samples. As a 
result, the approximation error is greatly reduced. The method according to the invention 
is low in computational complexity and modular in structure. 



25 
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What is claimed: 



1 . A method for transforming a digital signal from the time domain into the 
frequency domain using a transformation function comprising a transformation matrix, 
5 the digital signal is divided into a plurality of blocks, each block comprising a predefined 
number of data symbols, the method comprising: 

- transforming the two sub-blocks by a transforming element, the transforming element 
comprises a plurality of lifting stages, wherein the two sub-blocks are received by a first 

10 lifting stage of the plurality of lifting stages as first and second input signals, and wherein 
each lifting stage comprises a transformation path, wherein the first input signal received 
by the lifting stage is processed by a domain transformer and a rounding unit, and is 
subsequently summed with the second input signal in a summation unit to result in an 
output signal; and wherein each subsequent lifting stage receives the output signal of the 

15 previous lifting stage as the first input signal, and the first input signal of the previous 

lifting stage as the second input signal, to generate the corresponding output signal of the 
subsequent stage. 

2. The method of claim 1 , wherein the transformation function is a DCT-I 
20 transformation function, DCT-IV transformation function, DST-I transformation 
function, DST-IV transformation function, FFT-I transformation function, FFT-IV 
transformation function, DWT I transformation function or DWT IV transformation 
function. 

25 3. The method of claim 1 or 2, wherein each lifting stage corresponds to a lifting 
matrix, the lifting matrix comprises four sub-matrices with two permutation matrices as 
two of the sub-matrices in one diagonal, and with the transformation matrix and a zero as 
the omer two of the sub-matrices in the other diagonal. 

30 4. The-method of claim-3, wherein_the permutation matrices in each lifting matrix 
are identity matrices. 

V 

14 
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5 The method of any one of claims 1 to 4, wherein three lifting stages are used in 
the lifting ladder model for transforming the two blocks of the digital signal. 

6. A method for reversibly transforming a digital signal from the time domain into 
the frequency domain, the digital signal is divided into a plurality of blocks, each block 
comprising a predefined number of data symbols, the method, for a first and second block 
of the digital signals, comprising 



10 - generating an intermediate signal by 

-transforming the second block of the digital signal to a transformed second signal 

using a transformation function; 

- rounding the transformed second signal to a rounded second signal such that the 
value represented by the transformed second signal is rounded to an integer value; 

15 and 

- adding the first block of the digital signal to the rounded second signal, 

- generating a first output signal by 

- transforming the intermediate signal to a transformed intermediate signal using 

the transformation function; 
20 - rounding the transformed intermediate signal to a rounded intermediate signal, 

such that the value represented by the transformed intermediate signal is rounded 
to an integer value; and 

- subtracting the second block of the digital signal from the rounded intermediate 
signal^ 

25 - generating a second output signal by 

- transforming the first output signal to a transformed output signal using the 

transformation matrix; 

- rounding the transformed output signal to a rounded output signal, such that the 
value represented by the transformed output signal is rounded to an integer value; 

30 and 

- subtracting the rounded output signal from the intermediate signal, 
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wherein the first and second output signals are the transformed first and second block of 
the digital signal, respectively! 

7. A device for transforming a digital signal from the time domain into the 
5 frequency domain using a transformation function comprising a transformation matrix, 
the digital signal is divided into a plurality of blocks, each block comprising a predefined 
number of data symbols, the device comprising: 

- a transformation unit for transforming the two sub-blocks by a ti^sforming element, 
10 the transforming element comprises a plurality of lifting stages, wherein the two sub- 
blocks are received by a first lifting stage of the plurality of lifting stages as first and 
second input signals, and wherein each lifting stage comprises a transformation path, 
wherein the first input signal received by the lifting stage is processed by a domain 
transformer and a rounding unit, and is subsequently summed with the second input 
1 5 signal in a summation unit to result in an output signal; and wherein each subsequent 
lifting stage receives the output signal of the previous lifting stage as the first input 
signal, and the first input signal of the previous lifting stage as the second input signal, to 
generate the corresponding output signal of the subsequent stage. 

20 8. A computer readable medium, having a program recorded thereon, wherein the 

program is to make the computer execute a procedure transforming a digital signal from 
the time domain into the frequency domain using a transformation function comprising a 
transformation matrix, the digital signal is divided into a plurality of blocks, each block 
comprising a predefined number of data symbols, the method comprising: 

25 

. transforming the two sub-blocks by a transforming element, the transforming element 
comprises a plurality of lifting stages, wherein the two sub-blocks are received by a first 
lifting stage of the plurality of lifting stages as first and second input signals, and wherein 
each lifting stage comprises a transformation path, wherein the first input signal received 
30 Ky tfae lifting st age is processed by a domain transformer and a rounding unit, and is 
subsequently summed with the second input signal in a summation unit to result in an 
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output signal; and wherein each subsequent lifting stage receives the output signal of the 
previous lifting stage as the first input signal, and the first input signal of the previous 
lifting stage as the second input signal, to generate the corresponding output signal of the 
subsequent stage. 

5 

« 

9. A computer program element which is to make the computer execute a procedure for 
transforming a digital signal from the time domain into the frequency domain using a 
transformation function comprising a transformation matrix, the digital signal is divided 
into a plurality of blocks, each block comprising a predefined number of data symbols, 
10 the method comprising: 

- transforming the two sub-blocks by a transforming element, the transforming element 
comprises a plurality of lifting stages, wherein the two sub-blocks are received by a first 
lifting stage of the plurality of lifting stages as first and second input signals, and wherein 

15 each lifting stage comprises a transformation path, wherein the first input signal received 
by the lifting stage is processed by a domain transformer and a rounding unit, and is 
subsequently summed with the second input signal in a summation unit to result in an 
output signal; and wherein each subsequent lifting stage receives the output signal of the 
previous lifting stage as the first input signal, and the first input signal of the previous 

20 lifting stage as the second input signal, to generate the corresponding output signal of the 
subsequent stage. 
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Figure 2. Forward IntDCT-IV flowchart 
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Figure 3. Inverse IntDCT-IV flowchart 
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□ COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

□ LINES OR MARKS ON ORIGINAL DOCUMENT 

■ * 



UJ REFERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 
□ OTHER: 



IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 



BEST AVAILABLE IMAGES 



□ BLACK BORDERS 





